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Abstract 

The modular multiplication operator, a central subroutine in Shor's factoring algorithm, is shown 
to be a coherent superposition of two quantum bakers maps when the multiplier is 2. The classical 
limit of the maps being completely chaotic, it is shown that there exist perturbations that push 
the modular multiplication operator into regimes of generic quantum chaos with spectral fluctua- 
tions that are those of random matrices. For the initial state of relevance to Shor's algorithm we 
study fidelity decay due to phase and bit-flip errors in a single qubit and show exponential decay 
with shoulders at multiples or half- multiples of the order. A simple model is used to gain some 
understanding of this behavior. 
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I. INTRODUCTION 



Given an M dimensional complex Hilbert space consider an orthonormal basis |m), m = 
0, . . . ,M — 1. The modular multiplication by a number p coprime to M is the unitary 
permutation operator Up-. 

Up\m) — > \mp {mod M)). (1) 

Repeated application of Up is also the modular exponentiation operator and lead to Up\m) = 
|mp^ (modM)). Up is periodic, i.e. there exists ko such that f/^" = /, or = l(modM). 
This period though is an irregular function of M for a given p and is the "multiplicative 
order" ofp(modM) TMs operator . a cruciaUubroutine of Sho.s factonng algorithrn 
pl] in which it is required to perform the following operation on a bipartite Hilbert space: 

|j)|l)^|j)b^modM) (2) 

where < j < 2* — 1 and M < 2*. Typically M is a large number and therefore we 
need to calculate large powers and their residues modM. This is considerably simplified by 
modular arithmetic [3!] and the whole modular exponentiation step can be performed with 
(9((logM)^) number of gates Once this is done, a quantum Fourier transform extracts 
the period ko with reasonable rate of success. Given the period (order) it is possible to 
find a factor efficiently, provided that pk^ is even and p'^o/^ ^ — l(modM) by well-known 
procedures of number theory, and using only classical computers jsl. 

Shor's algorithm exploits the polynomial speed of the quantum Fourier transform to find 
the order and hence offers a polynomially scaling algorithm for factoring numbers. Shor's 
algorithm has been implemented experimentally |5| although the number so factored is still 
very small to excite any practical application. The effect of decoherence and gate errors on 
Shor's algorithm are irnportant considerations and have been addressed by several authors 



previously 
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lit]. For instance Ref. ^discuss the impact of environmental 



decoherence on the algorithm, while in [9(] direct detailed simulations have shown that the 
Shor algorithm is highly sensitive to gate errors, and the effect of static imperfections have 
been studied recently in[l]J]. In this paper we will not directly simulate Shor's algorithm but 
look closely at the modular multiplication and exponentiation for the special and simplest 
case p = 2. We will call U2 as S, the shift operator as it performs the simple action of a 
qubit cychc shift if M is a power of 2. 
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Quantum algorithms have been studied earher with a view to see if they had properties of 
quantum chaotic systems [l^. Recently it was shown that the spectrum of the unitary 
jart of Shor's algorithm, properly desymmetrized had typical random matrix fluctuations 



131], indicating that the operator itself may be quantum chaotic. However it was also 



pointed out that the origin of the chaos is the modular exponentiation part which is akin 
to nongeneric quantum chaotic systems such as the cat maps [l4^. Here we make this 
connection more precise and show that the classical limit of these subroutines is an ad- 
mixture of two bakers maps. Bakers maps are paradigms of deterministic classical chaos 



that random as a coin toss 
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17l | . The dimensionless inverse Planck constant in the 



Shor algorithm is the number to be factored and hence the classical limit is reached through 
a practically important regime. We show that due to the proximity with such operators, 
there are perturbations that push the modular exponentiation part (and therefore indeed 
the whole of Shor's algorithm) into regimes of generic quantum chaos. However we are 
not that much interested in stationary state properties as in time-evolving states, in fact 
on those states that are used in Shor's algorithm. Hence we study the fidelity of repeated 
modular multiplication, or the modular exponentiation, and show how the decay depends 
on the classical limit. We also provide a simple model for the exponential fidelity decay that 
is exactly solvable and captures the actual behavior reasonably well. 



II. THE BAKER AND THE SHIFT OPERATOR 



The classical baker's map 



16|, 



17( 1 Be is an area-preserving transformation of the unit 



phase-space square [0, 1) x [0, 1) onto itself, which takes a phase-space point {q,p) to {q',p') 
where {q' = 2q, p' =p/2) if < g < 1/2 and {q' = 2g-l, p' = (p+l)/2) if 1/2 < g < 1. The 
stretching along the horizontal q direction by a factor of two is compensated exactly by a 
compression in the vertical p direction. This is well known to be a fully chaotic system that in 
a mathematically precise sense is as random as a coin toss [18]. The area-preserving property 
makes this map a model of chaotic two-degree of freedom Hamiltonian systems. The lack of 
a generating Hamiltonian is compensated by the existence of a classical generating function 
of the canonical transformation Be- The chaos is inferred by expressing a phase space point 
in the binary representation, if g = 0.aoaia2 ■ ■ ■ and p = 0.a_ia_2a-3 ■ ■ ■ , where are either 
or 1, then q' = 0.aia2 ■ ■ ■ and p' = 0.aoa--ia-2'i-3 ■ ■ ■ • Thus the most significant bits 
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of q are lost at the rate of a bit per iteration, leading to an exponential increase in any 
initial error. The Lyapunov exponent is log(2) per iteration. This "left-shift" is in fact an 
important mechanism for the generation of Hamiltonian chaos 17|] , and in more complicated 
forms arises generically. 

This was quantized first by Balazs and Voros[l^, 20|. Quantization in this context is the 
construction of an appropriate unitary operator that evolves states over one iteration and has 
the correct classical limit. Symmetries that may be broken on quantization must be restored 
in this limit. The Hilbert space of states is finite dimensional, has position and momentum 
states, denoted by |g„) and \pm). If periodic boundary conditions are assumed, \qn+N) = 
\pm+N) = \Pm), this implies that the transformation functions between position and 
momentum is the discrete Fourier transform: {Fi^)mn = {Pm\<ln) = exp[— 27rzmn/A^]/v/iV, 
m, = 0, 1, 2, . . . , — 1. Here is an effective scaled Planck constant as N = A/h = 1/h, 
where A is the area of the phase space, here unity. Thus the classical limit is the large 
limit. If B is the quantum baker's map, Balazs and Voros required that {pm\B\qn) = 
{Pm\q2n) = {FN/2)mn if ^ and 111 are both < N/2 — 1. This is almost like requiring that 
B takes to \q2n) mimicking the classical stretching action, except that the momentum 
components above A^/2 are set zero. {{pm\B\qn) = for p^ > A^/2 and g„ < N/2). It is also 



clear from this that B is very close to the action of modular multiplication with p = 2 |21 |. 
A^ is throughout assumed to be an even integer. In fact we will set N = 2L and can then 
consider the quantum baker to act on a Hilbert space of a qubit coupled to an L dimensional 
systems. 

A similar argument is made for the second half of the transformation, and remarkably 
these conditions are consistent and produce an unitary operator which has a broken parity 
symmetryc |2o|. The classical symmetry being {q — > 1 — q, p — > 1 — p). Saraceno 22] 
restored this by imposing anti-periodic boundary conditions, and this leads to the quantum 
baker's map: 

^ f Gl \ 

B = G,l[ ' ], (3) 
\ GlJ 

where {GN)mn = {Pm\qn) = exp[— 27ri(m + l/2)(n + 1/2)/N]/^/N. B is an unitary matrix, 
whose repeated application is the quantum version of the full left-shift of classical chaos. This 
quantum map has been continued to be studied as it has many properties of generic quantum 



chaotic systems, including random matrix like spectral fluctuations 



20| and eigenfunction 
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scarring [22|]. It is also amenable to a simple semiclassical periodic orbit sum, and hence 



has been used in the study of such approximations 
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2J]. For that are powers of 2 



it was found that the Hadamard and related transforms highly simplified the eigenstates 



and some of them are remarkably well described by the Thue-Morse sequence 
its Fourier transform 
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26| and 



27l |. It has also been used in the study of entanglement [28 1 



and hypersensitivity of quantum chaos [29|. It is possible to design a quantum circuit for 



the quantum baker map [s^ and this been implemented on a NMR quantum computer 
experimentally |3l| . 

If one is not mindful of classical symmetries being fully preserved on quantization, there 
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32| . An important class of such 



are a large number of possible quantum baker maps 
"decorated bakers" are got by embellishing the original bakers map with relative phases 
in the half-sized Fourier blocks, as well as in the definition of the Fourier transform itself, 
as done below. In a previous work we have constructed such a decorated quantum bakers 
map using the shift operator S 2l|]. It will be useful to do the converse and construct the 
shift operator from the quantum bakers map or similar operators. It is well-known that 
the "hard" part of Shor's algorithm is the implementation of the modular exponentiation 
step. On the other hand the quantum bakers map is implemented with quantum Fourier 
transforms (QFTs) and this may make the implementation of the shift operator possible with 
the QFTs. We explicitly show this at least for the case p = 2. More importantly for us it 
will enable embedding the shift operator in a larger family of operators which includes maps 
with well defined classical limits, thereby making the classical limit of modular multiplication 
explicit. 



The shift operator we have already defined, however we restate it for clarity as: 



S\n) = |2n(modiV - 1)) 



(4) 



with the caveat that 5*1 — 1) = |A^ — 1). This corresponds to our earlier definition with 
M = A^ — 1 with one more state (|A^ — 1)) added to the Hilbert space, but which remains 
fixed, and outside any dynamics we are interested in, but may participate when there are 
perturbations. Note that since for the bakers map A^ is an even integer S is unitary. We 
define a generalized Fourier transform as 

1 



exp(-27r(n + a)(m + /3)/N) 



(5) 
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Evaluating the product of the Fourier transform and 5* we can derive, merely by summing 
finite geometric series and using elementary properties of exponentials that 

S = ^F^L (a, ") , I • (6) 




Note that the operator S does not depend on the phase a that appears on the R.H.S.. We 
use this freedom to break or keep the parity symmetry. A natural and simple choice is a = 0, 
but a = 1/2 leads to symmetric operators as explained below. The structure of the above 
identity allows this to be written as 

+ 5k) (7) 



where 



V2 



B,, = F-,\a,a) | \' . | . (8) 



and 



B',, = F-,\a,a) | . ^ '^"^ |. (9) 

e-*""FL(a,f) 

Thus remarkably the modular multiplication S can be written as a sum of two unitary 
operators (with a normalization factor). B2L and are two quantum baker maps, of which 



the former one is the standard one, which we have discussed above. It is well 
such decorated bakers also perform the same classical actions as normal bakers 



inown that 
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32], the 



classical limit being L — > 00. The operator i?2i has not been studied nearly as much. 



but has recently appeared in a work that uses this to study coupled chaotic systems [34 1- 
The classical limit (say B'^) as pointed out in this study corresponds to a different stacking 
order of the vertical partitions of the bakers map after they have been stretched. Instead 
of the usual left-half transiting to the bottom-half, it is put in the top-half and the right- 
half goes into the bottom half. This fixes the lower right-hand corner of the square. Thus 
B'^{q,p) = (2g, {p + l)/2) if g < 1/2 and (2g - l,p/2) if g > 1/2. Again the operator B'^^ 
that appears above differs from the one used earlier in terms of the "decorations". Of course 
these decorations are absolutely crucial so that the two unitary evolutions, which are non- 
periodic and have random matrix like properties, add and conspire to produce the simple 
shift operator that is completely periodic. Previous studies of the classical limit of operators 



such as S include those of what is called the "extremal quantum baker map'' 



35] and it has 



FIG. 1: A schematic view of the shift operator as a sum of two bakers maps, one with the usual 
stacking order of the vertical left half being stretched to the bottom horizontal half (i?c-classical, 
-B2L-quantum) and the baker with a reverse stacking order (i?^-classical, i32/^-quantum) . 



been suggested that the classical limit corresponds to a "stochastic classical map" [36|, |3j 
In Fig. ([1]) we have shown a schematic of the classical bakers maps that on quantization and 
coherent addition yield the shift operator. Also see Ref. [21|] for a description and figure of 
the action of S on Weyl coherent states. 

That the simple shift operator's can be thought of as a coherent superposition of two 
quantum chaotic evolutions has been demonstrated above in a particularly simple way. This 
suggests that there maybe perturbations of the operator S that ar e g eneric and may possess 



random matrix [38] like and other quantum chaotic properties 39|, |40]. We show below 



that this is indeed the case. Since quantum chaotic operators also are typically sensitive to 



perturbations 



42|, this may have implications for the operation of the Shor algorithm. 
We partly study this by measuring the fidelity of 5* to small perturbations and show that 
the fidelity decays exponentially in time till the order of — 1 or half of this. Thereafter 
it typically shows enhanced rate of decays at multiples of this time, but could also show 
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strong recurrences. Surprisingly a simple analysis when N is a, power of 2 captures many of 
the qualitative features of the more general case. 



III. PERTURBATIONS OF THE SHIFT OPERATOR AND QUANTUM CHAOS 

In terms of operations on the Hilbert spaces of the tensor product 'H2®'Hl we may write 
the shift operator as 

1 / 1 I \ ( FAa,^) \ 

S = F-l{a,a)o—\ \®Ilo[ ' . (10) 

Thus the modular exponentiation maybe implemented with QFTs that have suitable phases. 
However the dimensionalities of the QFTs are not in general powers of 2 and are therefore not 
the standard ones in use. We choose to perturb the central operator in the above equation, 
and perturb only the qubit space 7^2- In particular we consider the smooth embedding of 
the shift operator in the family: 

1 / 1 I \ ( FAa,^) 

V2 \^e-*'^" -e-'""^ J y Fc(q;, i±^) 

= F-j^(a, a) o exp(-i^P) ® 7^ o F2l(q;, a) S ^ V(e)S. (11) 

here P is the perturbing Hermitian operator on the qubit space, and V{9) defined through 
the last equation is the Fourier transform of the perturbation generated by P. The operator 
5'(0; |, P) is the unperturbed shift operator simply called S so far. The family of operators 
S{9;a,P) now depends on the phase a as well, although S{0;a,P) does not. We display 
this dependence explicitly as the phase a does play a crucial role. 

The operator S has the quantum parity symmetry R : R\n) = \N — n — 1) , that is 
SR = RS. If L is a power of 2, then R is simply the product ^'^^a^. Perturbations 
of the shift will therefore in general approximately preserve this symmetry. To analyse 
random matrix properties it is desirable to completely break a symmetry or preserve it and 
desymmetrize the operators. Since we want to retain the character of a small perturbation, 
we first choose to preserve the parity symmetry exactly. We can do this by adopting anti- 
periodic boundary conditions, a = 1/2 and choosing P = ax- This will lead to the family: 

S(e-\) = F~'(^ ^)( ^i"(f-^)^^(i'i) cos(f-^)F4i,|) 
2' -^2'2M _,eos(f-^)F,(i,i) i sin(f - ^)F,(i |) 



When 9 = ±7r/4, the operators correspond to the bakers of type and B2L respectively. 
For other angles it represents a coherent mixture of the two types of baker operator stacking 
while at ^ = it is the usual shift operator. 

We take the even subspace of the spectrum of S{6; |, a^) and show in Fig. ([2]) the nearest 
neighbor spacing statistics for two case of small angles 6. It is clear that if the perturbation 
is very small, the rigid, harmonic oscillator like spectrum of S widens into one where there is 
dominant level repulsion, and for fairly significant perturbations becomes a generic one that 
belongs to the universality class of the Gaussian Orthogonal Ensemble (GOE) of random 
matrices 38 L well-known to apply to quantum chaotic systems that have time- reversal 



symmetry 40|, |43|. Previously it was shown that the unitary part of the full Shor algorithm. 



inclu ding the Hadamard and the Fourier transforms, had fluctuations that were of the CUE 
kind '13]]. We note here that restricting ourselves to the modular multiplication part with 
a particular perturbation allows us to preserve the time-reversal that holds for individual 
quantum bakers maps [20|]. 

Thus we see that indeed there are perturbations of the shift-operator that are quantum 
chaotic. There is also a crucial dependence on the number (or the number to be factored 

— 1). If were a power of 2 such as 4096, instead of 4094 in Fig. ([2]) there would 
be much more deviation from the GOE distribution, with a large peak near the origin. 
This anomalous statistic arises from the extreme degeneracy of the eigenangles when is 
a power of 2, and is special. A similar situation where there is an extreme dependence 
of the statistics of the spectrum on the effective Planck constant A^ arises in the case of 
the perturbed cat maps 4J], and presumably for similar reasons. Earlier it has also been 
pointed out that perturbing the cat maps slightly so that the sawtooth map arises leads to a 
rapid restoration of the generic fluctuation characteristics of quantum chaotic systems j45|. 
Thus the similarities of the shift map to the quantum chaotic cat maps with their special 
arithmetic properties is further highlighted here. 



IV. FIDELITY DECAY 



We turn to non-stationary properties, as indeed the Shor algorithm is the result of time 
evolution of a particular initial state which corresponds to the state The algorithm 
requires requires finding the states mod (A^ — 1)) for x coprime to A^ — 1. As stated 
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0.5 1 1.5 2 2.5 3 
Spacing s 

FIG. 2: The nearest neighbor spacing distribution of even-subspace eigenangles for two perturba- 
tions of the shift operator when A'^ = 4094, the perturbations preserving the parity symmetry. The 
smooth curve shows the corresponding GOE result of random matrix theory, the Wigner-Dyson 
distribution. 

previously we take x — 2 throughtout, and we now study how gate errors would proliferate 
in time. In particular we study the fidelity 

i{t) = mS-'S\e-a,P)\\)\^ = |(2*mod(iV-l)|5*(^;a,P)|l)|2 (13) 

In this section choose o; = for simplicity and note that there is a weakly broken parity 
symmetry as a result of this. 



A. Case: P = Ux 

The first case we take will be a rather special one wherein the perturbation is ineffective: 
the bit-flip P — a^- In this case we have that 

^e--F,(0,0) >-F,(0,i)\ 
^e-^F40,0) ^e^^F40,i) ) 

and f{t) — 1 for all time t. Note that there are only phases multiplying the Fourier blocks, 
and therefore the classical hmit of this family of operators is the same as that of the simple 
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5(^;0,a,)=F-/(0,0) 



shift operator: a coherent sum of two baker maps. However, the reason the fidehty is unity 
for all time is due to a rather intriguing if simply verifiable identity. We will write F2L for 
-^2l(0,0) below. The perturbation operator is: 

y{d) = F^L (exp(-20(T,) h) F2L = h® exp(-zea,). (15) 

The last equality is an identity valid for all integer L. This in turn simply follows from the 
identity: 

F^l (a, ® II) F2L = {II ® a,) (16) 

which maybe directly verified. Note that since cr^ ® /l is a circulant matrix it has to be 
diagonalized by the Fourier transform, and since the eigenvalues are ±1 these are the only 
possible diagonal entries. It is also easily verified that this has the structure of L repetitions 
of (1,-1) pairs which are the diagonal entries of a^. In some sense the Fourier transform 
is simultaneously performing a bit reversal and a ninety degree rotation in qubit space. 
However note that for this identity to be true we do not require that L be a power of 2. We 
state here that similar identities do not hold for the other two Pauli matrices, but there are 
approximations that we will state further ahead. Since az merely changes the phase of the 
state (in standard basis) it follows that the fidelity f(t) = 1 always. 

B. Case: P = ay 

In this case we see that the classical limit is altered by the perturbation. The Hadamard 
transform in the qubit space is further rotated around the y-axis in spin space and the final 
operator is similar to that used in Eq. ( |T2l) . which we recall is for the case when the phase 
a = 1/2 and for a (7^ perturbation. 

\^cos(J-»)fi(0,0) -sm(J-»)Fi(0,i) j 

Therefore this case is the closest to the parity preserving case we have already discussed 
and shown the sharp transitions to features of a quantum chaotic spectrum. In Fig. ([3]) we 
plot the fidelity /(t) for a set of values that are close to 256. We notice immediately 
that although the N values are as close as can be (A^ must be even) the fidelity decays in 
qualitatively different manners. Except for a very short-time scale the decays are different 
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I ^ ^ 1 — -"'"^ 

50 100 150 200 250 

time t 

FIG. 3: The fidelity decay for four neighboring values of N. The perturbation is P = ay, the phase 
a = 0, and 61 = 0.05. 

and one sees a prominent "shoulder" in each of the curves at which the fidelity starts to 
decay even faster. 

The easiest case to discern this in the figure is for N — 254 when the shoulder occurs 
at t — 110. It is quite easy to numerically relate the time at which this occurs to the 
multiphcative order of 2 modulo N — 1, referred to henceforth loosely as simply the order of 
N-1. Recall that this is the smallest number A;o(A^-l) such that 2^o(^-^) = 1 modulo(iV-l). 
Wc arc guaranteed that such a number exists because — 1 is an odd integer. Indeed 
A;o(253) = 110, while A;o(255) = 8, A;o(251) = 50, and A;o(249) = 82. Thus the fidelity shows 
a shoulder either exactly a,t t — ko{N — 1) or a,t t — ko{N — l)/2, the first case is observed 
for N — 256 and N — 254 while the latter is the case for N — 252 and 250. It is significant 
then that for numbers of larger orders the fidelity can decay considerably even for small 
perturbations. Note that 9 = 0.05 in the figure which roughly translates to a 0.51 : 0.49 
mixture of the two types of bakers, while a 0.5 : 0.5 "mixture" will be the unperturbed shift 
operator. The larger the order is the higher powers of S must be calculated and the higher 
chance of the fidelity to be lowered. It is interesting that the objective of the Shor algorithm 
namely finding the order already appears in the fidelity as a critical time. 

We can gain a qualitative understanding of these behaviors with a surprisingly simple 
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model. Consider the case when is a power of 2 say = 2 and let the perturbation be 

V{9) = Il® exp(-i^a^). (18) 

Note that this has the same structure as the perturbation from the previous case, but is not 
the true perturbation in this one. Then the initial state is |0 ■ ■ ■ 01) and it is clear that for 
t < M: 

{V{e)S)^ |0-- - 01) = l2M-t ® exp(-i^a^) ®exp(-i^a^) ■■ - ^ exp(-i0a^)|O- - -01). (19) 

and f{t) = |cos^(^)|*. Thus the initial fidelity decay is in fact exponential with a rate 
— log(| cos^(0)|). However beyond t = M there is an additional error adding up and so 
f{t) = \ cos2(^)|2^^-*| cos2(2e)|*-^^ for M < t < 2M. Thus in this range of time the fidely 
decays about four times as fast. We can write in general for this model that 

fit) = I cos2(r0)|(^+^)*^'*| cos2((r + l)^)]*"^^ (20) 

where r = [t/M] and [x] is the integer part of x. We have shown in Fig. (jlj) how good an 
approximation this can be for the case 2 situation when is indeed a power of 2 (or the 
number we want to factor is one less than a power of 2). There is even good quantitative 
agreement. On the other hand when N = 254, (253 has a high order of 110) the approximate 
formula is only qualitatively correct as seen in Fig. (jl]). 
The model works reasonably well because 

F-j^ (expi-iOay) ® h) F2L ^ (h ® exp(-i0a^)) , (21) 

which follows from 

{ay ® h) F2L ~ ih ® CT,) , (22) 

so that exact perturbation operator which would have been the L.H.S. of Eq. fl2Tl) is ap- 
proximated by its R.H.S. which we have used above. This is the counterpart of Eq. (fT6ll . 
however here this is only an approximation. 

This "model" or approximation does not explain the appearance of half of the periods 
for some values of A^, such as for 252 and 250 above. Indeed when A^ is a power of 2 we 
will always observe the first shoulder at the order of A^ — 1. This is in fact the result of 
the possibility that there exists an integer k' such that 2^' mod (A^ — 1) = —1, which implies 
that k' = ko/2. If there exists such an integer then we must, according to Shor's algorithm, 
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N=256. Exact values □ 
Approximate model 




40 60 
time t 



N=254. Exact values □ 
Approximate model 




200 



FIG. 4: The comparison of the fidehty decay for the ay perturbation (a = 0, and 9 
the analytical estimate in Eq. pop from an approximate model for two cases of N 



0.05) with 



choose a different integer (other than 2) to find its order of. That is we cannot use the order 
of 2 to find a factor of A^ — 1, which is the ultimate objective. For our analysis, this situation 
implies that 



I 



\ 



1 



R 



N~2 



\ 



V 



(23) 



where Rn-2 is the parity operator with 1 along its secondary diagonal and zero elsewhere. 
That is 5'''=o(Af-i)/2 jg ajj^ost the parity operator except that and N — 1 instead of being 
interchanged are fixed by S, and hence all its powers. To clarify Eq. fl23l) may or may not hold 
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depending on if 2'^o(^~^)/^ = —1 mod (A^ — 1) or not. For instance this is never the case if 
is a power of 2. We note that this simphfied model has also been considered by 29| recently 
to show that exponential fidelity decay does not necessarily mean an hypersensitivity to 
perturbations. However in the context of this paper it is interesting that the model works 
approximately even when is not a power of 2 and preliminary results indicate that there 
is hypersensitivity to perturbations as well 46|. 
In general for an initial state {ipo) "we have that 

/(t) = |(^o|V^tV^*-i---^i|^o)r (24) 

where Vi is the perturbation in the interaction picture: VJ = S~''VSK Thus always Vkg(N-i) = 
V but in the situation where Eq. ( 123|) holds we have that Vko{N -i)/2 = R'~^VR'. Thus it 
is clear that these times are special for the fidelity as seen in the numerical calculations as 
well. While these arguments along with the approximate model gives a fair understanding 
of the decay, it is not complete and a more detailed analysis of the product above must 
be carried out, which the author is unable to provide. There is a rather large literature 
surrounding the so-called Loschmidt echo [47] , or fidelity, wherein quantum chaotic systems 
have been subjected to a small perturbation on reversal. The current discussion is in fact 
closely related, but previous work has naturally concentrated on the generic case of a non- 
degenerate operator that is perturbed. In the case of the shift operator, it can be highly 
degenerate, as well as completely periodic, thereby making it "non-generic". It can be 
compared again to the quantum chaotic cat maps that are also periodic and degenerate in 
general. Smooth perturbations of this for instance of the type that has been studied before 
50l | could produce fidelity decays of a similar character that we have noted here. 



C. Case: P = 

In this case 



^e^^Fi(0,0) ^e^^Fz.(0,i; 



(25) 



2' 
2' 

Note that there seems to be only a minor change, namely those of signs of phases in the 
Fourier blocks, compared to the first case, and also the classical limit still remains unaltered 
by the phase-flip perturbation. However, the fidelity does decay even due to the "quantum 
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FIG. 5: The fidelity decay for four neighboring values of A^. The perturbation is P = cj^, the phase 
a = 0, and 6* = 0.05. 

perturbation" and seems to be of a similar character to that observed when P = ay, namely 
the previous case. The differences start showing up sharply in the case when Eq. f l2^ holds, 
namely when ko{N — 1) is such that 2''°^^'^^ = —1 mod (A^ — 1). It appears to be generically 
the case that beyond this time there are large oscillations reminiscent of fidelity decay in 
near-integrable systems 4^. This is illustrated in Fig. ([5]) for values of A^. In the cases 
when = 252 and 250 Eq. (1231) holds and we see that beyond time of half the order there 
are regular oscillations with this period. 

That we must expect a fidelity decay is due to a counterpart of the approximation used 
in the previous case, namely 



(26) 



which follows from 

F^-/ (a, ® h) F2L ^ ih ® cxy) , (27) 

which is a result of combining the identity in Eq. ( |T6l) and the approximation in Eq. ( l22l) . 
Thus this final case of perturbation we consider is sort of intermediate between cases A and 
B, however for practical purposes it is closer to case B, as the time behavior beyond the 
time of the order or half the order is not likely to be of interest from the point of view of 
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the Shor algorithm. The approximate formula in Eq. (120|) continues to be approximately 
good for those for which the condition in Eq. fl23l) does not hold, and for those for which 
it does, it is approximately good till time of half the order at which the oscillations begin. 



V. DISCUSSION 

We have studied three archetypal perturbations, phase-flip, bit-flip and a combination 
therefore, that are possible in the critical part of Shor's algorithm, namely the modular mul- 
tiplication or exponentiation part. We have confined ourselves to the simplest possible case 
when the multiplier is 2, when these perturbations can be interpreted in terms of coherent 
superpositions of quantum bakers maps, whose classical limits are completely chaotic, and 
are models of randomness. Thus we have shown that there are generic perturbations of the 
modular exponentiation operator that will qualify as "quantum chaotic". We have shown 
this by computing the nearest-neighbor spacing statistics and seeing that it is of the type 
expected of random matrices. More pertinent to the algorithm itself we have studied the 
fidelity decay that occurs with the relevant initial state and shown that for the three types 
of perturbations there are three possible fidelity decay behaviors. This can be interpreted in 
terms of the fact that some perturbations alter the classical limit while some do not, as well 
as in arising from some identities (one exact and one approximate) that involve the Pauli 
spin matrices and the Fourier transform, which while the author has not seen before, are 
completely elementary and likely to be known and useful already. A simple model of the 
fidelity decay is afforded by these identities that describes surprisingly well the exponential 
decay in time punctuated by shoulders at times related to the order. An exact solution 
of the problem seems unlikely, and semiclassical analysis cumbersome due to the fact that 
the modular exponentiation (when the multiplier is 2) is essentially the sum of two unitary 
operators with well defined classical limits. 

The precise impact of the exponential fidelity decay on the functioning of the algorithm 
remains to be seen. Such an study for the case of static imperfections was recently carried 



out ll|. We have been primarily interested in pointing to the deep and exact relationship 
between the modular exponentiation part of the Shor algorithm and the quantization of an 
archetypal model of classical deterministic chaos, namely the bakers map. If we had larger 
multipliers that 2, as will generally be the case, it is reasonable to expect that these will 
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be related to generalized bakers maps with more than 2 partitions [20|]. The number of 
possible stacking are also more, but it is completely conceivable that once again there are 
perturbations to the modular exponentiation operator that are quantum chaotic and close to 
the quantization of such bakers. That such bakers will have larger Lyapunov exponents and 
have greater classical randomness may make the quantum operators even more susceptible 
to such gate errors; however this is at the moment mere speculation. 
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